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Abstract 

In this manuscript we analyse the behaviour of the probability density function of the sum 
of N deterministic variables generated from the triangle map of Casati-Prosen. For the 
case in which the map is both ergodic and mixing the resulting probability density function 
quickly concurs with the Normal distribution. When these properties are modified the 
resulting probability density functions are described by power-laws. Moreover, contrarily 
to what it would be expected, as the number of added variables N increases the distance 
to Gaussian distribution increases. This behaviour goes against standard central limit 
theorem. By extrapolation of our finite size results we preview that in the limit of iV 
going to infinity the distribution has the same asymptotic decay as a Lorenztian (or a 
q = 2- Gaussian). 



1 Introduction 

The central limit theory has been a subject of studied within the natural sciences for many 
generations. We might even state that CLT originated in 1713 with Bernoulli's weak law 
of large numbers |lj. After Bernoulli, de Moivre [2], Laplace [3], and Gauss, amongst 
others, made crucial contributions to the establishment of the Normal probability density 
function (PDF) as the stable distribution when one considers the sum of independent and 
identically distributed random variables with finite standard deviation. The stability of 
the Normal distribution, the central limit theorem (CLT), was just formally established 
by the Russian mathematician Lyapunov 188 years after Bernoulli [I]. Afterwards, Levy 
and Gnedenko E] generalised the CLT to account for independent and identically 
distributed random variables but with infinite standard deviations, followed by broader 
generalisations which include dependency between variables [TJ [S]. With the advent of 
computation in the 1970s, chaos theory and non-linear phenomena achieved huge progress. 
It was then possible to verify the existence of a central limit theory for deterministic 
variables as well [9| [TP]. More recently, central limit theory has been the focus of renewed 
interest within statistical mechanics mostly because of the endeavour to establish the 
optimising PDF of non-additive entropy S q [UJ as the stable distribution for the sum of 
random or deterministic variables for some special kind of correlation [12] , on the edge of 
chaos [TBI [TJ], or even for systems in a metastable state [Hj] . 

In the sequel of this manuscript we communicate results on numerical investigations of 
the distribution of deterministic variables which arise from the sum of variables generated 

^mail address: Silvio.Queiros@unilever.com, sdqueiro@googlemail.com 



1 



from the triangle map introduced by Casati and Prosen [16J. Our analysis is performed 
at two different regimes: a first illustrative one, for which the system is both ergodic and 
mixing, and a second in which the system is weakly ergodic and apparently weakly mixing 
H. In the former case, the convergence towards the Gaussian distribution is clear and 
easily explained according to the standard CLT. In the latter case, we detect an anoma- 
lous behaviour characterised by non-Gaussian distributions which become more distant 
from the Gaussian distribution as the number of variables added increases. If we take 
into account the behaviour shown for phase occupancy observables in other conservative 
dynamical systems, this result is against all odds and it emphasises the very particular 
properties of this dynamical system. The remaining of the manuscript is organised as fol- 
lows: in Sec. |2] we introduce the triangular map and some of its properties, and in Sec. [3] 
we present our numerical results. Some conclusions and remarks are made in Sec. HI 

2 The triangle map 

The triangle map, z t+1 = T (z n ) introduced by Casati and Prosen [T7J OS], corresponds 
to a discrete transformation on a torus z = (x,y) G [—1, 1) x [—1, 1) with symmetrical 
coordinates, 



where sign a; = ±1, and a and (3 are real parameters of the map. Function (mod* 2) 
represents a modified definition of mod 2 in the interval between —1 and +1. This map 
has emerged from studies on the compatibility between linear dynamical instability and 
the exponential decay of Poincare recurrences. Map (0Q) is parabolic and area preserving. 
For the Jacobian matrix we have, det J = 1, and its trace, Tr J = 2. 

Concerning the relevance of parameters a and /3, it is known that when both of the 
parameters are irrational numbers, the map is ergodic [16J. Moreover, it attains the 
ergodicity property, i.e., averages over time equal averages over samples, very rapidly. 
For such a set of values the map is also mixing [T7j in the sense that it has a continuous 
spectral density besides the property previously defined. Evaluating Poincare recurrences, 
it has been found that the probability of an orbit to stay outside a specific subset of the 
torus for a time longer than t behaves as exp [—/j, t], with /i being very close to the Lesbegue 
measure of that subset. This fact is in accordance with a completely stochastic dynamics. 
The exponential decay leads to a linear separation of close orbits which has been related 
to non-extensive statistical mechanics formalism via a generalised Pesin-like identity that 
bridges a g-generalisation of Kolmogorov- Sinai entropy [TU] and a g-Lyapunov coefficient 
from the sensitivity to initial conditions [20] • The entropic index of this map was found 



2 By strongly ergodic we mean systems for which the exploration rate r (t) = n (t) /N [{n (t) is the 
number of cells visited by a orbit until a discrete time t and N is the total number of cells (N — > oo)], 
averaged over several initial conditions, presents the same functional form as the exploration rate of a 
random model r^M (t) — 1 — exp [—t/N]. Concomitantly, by strongly mixing we refer to systems whose 
correlation function, C (t) decays much faster than with fi > 1. Weakly ergodic/mixing system do 
not follow each of corresponding the features [TH] ■ 
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When a = and j3 is irrational, the map is still ergodic [22], but is never mixing 
[23] • For the case (3 = two situations might occur [TU]. If a is a rational number, then 
the dynamics are pseudo-integrable and confined, whereas for a irrational, the dynamics 
are found to be weakly ergodic, with the number of y n taken by a single orbit increasing 
as InT (0 < n < T). Furthermore, the ultra-slow apparent decay of the correlation 
function measured for this condition does not provide sufficient evidence of mixing. In 
recent years, some analytical attempts in order to characterise the map Eq. ([T]) have been 
made [21]. Despite this, it has not been possible to prove or reject the mixing property 
for the Casati-Prosen map although strong indication of a non-mixing property prevails. 



3 Results 

In this section, we present the numerical results of our study. Namely, we have considered 
variables Xjv and Yn that are obtained from the addition of x and y variables of map (JTJ), 

N 



X N = ^ Xi , (2) 



i=i 

N 



Y N = Y,Vi> ( 3 ) 

i=l 

that we analyse after detrending and rescaling, u' N = (un — (un)) /&n (u renders either 
X or Y and is the standard deviation, a N = ^ X)i=i u i ) ■ 111 our surve Y we have 
neglected the case a = since it destroys the dependence of y on x and we have focussed 
on the following situations, 

case l J «= I K (V5- 1) -e- 1 ] , 1 

which corresponds to the ergodic and mixing case studied in Refs. [T6l [21] and 

case II : ja = rc~^, (3 = o| , (5) 

where the map is weakly ergodic. For each analysis, we have randomly placed a set 
of initial conditions I (typically 10 7 elements) within the torus z and we have let the 
map run. The probability density functions P (Xn) and P (Yn) are then obtained from 
these I initial conditions. Our numerical calculations have been performed using the 
Mathematica™ kernel which assures a symbolic computational procedure. Although 
analytical considerations in case I are in principle possible, we have skipped them because 
we have used case I as a benchmark of the peculiar behaviour of case II that we are going 
to present. It is also worth stating that, for case I, where analytical considerations have 
been made (or can be made), conditions of strong mixing and (semi) conjugation to a 
Bernoulli shift are mandatory [H]. These conditions are not verified in our case II. 
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For case I, as expected from ergodic and mixing properties of map (pQ), both X N and 
Y N (detrended) approach the Normal distribution [TP] , 



G ( Un ^ = ^ — r exp 

V 27TCT N 



(6) 



as N goes to infinity, see Fig. CD Moreover, as it occurs for the Lyapunov CLT, cxtv follows 
the scaling relation, 

a N = y/Na u (7) 

with <Ti = 4^ as in Fig. [21 In Ref. [TU] a similar kind of analysis has been made by 
considering a cylinder y G (—00,00) 

y n = Vo + (3n + ap n , (8) 

with p„, G Z. 

We have also verified a skew in our PDFs, for small N, which are not visible in the 
PDFs of Ref. [IH], but might be comprehended according to analytical work made on 
other types of maps [HI EE] ■ Explicitly, skewed distributions have analytically been found 
when studying the same problem using the dissipative, fully chaotic and strongly mixing 
map Xt+i = 1 — 2 x\ [13] . 

A completely different behaviour is found when case II is analysed. For this case, 
we have concentrated on X^, although for Y/y we have obtained the same qualitative 
results. Instead of distributions reminiscent of a Normal distribution, we have numerically 
observed probability distributions which are well described by, 

P{u) ~ M-"- 1 (l«!i<iV), (9) 

with rj < 2 for every value of iV analysed, see Tab. [1] and Figs. [6] and HI The method 
applied to determine 77 has been the Meerschart-Scheffler estimator (see Appendix) [25] . 
The subsequent application of the Hill estimator [26] has given concordant results. 

The upper bound of 77 we have found (77 = 1.80) imposes that the standard deviation 
would diverge if the variable X(Y) was defined over the whole interval of real numbers. 
Since we are treating cases for which iV is finite, the support of the resulting PDFs is 
compact and defined between —N and iV for Xjy and Yn- This obviously leads to a 
finite standard deviation, 0^, for the cases we have studied. We have verified that <jn 
does not follow the scaling relation Eq. j7j), see Fig. [SJ Instead, a power-law dependence 
can approximately be given with an exponent close to 0.94 ± 0.02. In addition, we have 
observed that the shape of the distribution P (Xn) has strong similarity with the a-stable 
Levy distribution 

j p+00 

C a (Xn) = — / exp [— i k Xn — a \k\ a ] dk, (10) 
2 vr ./-oo 

((0 < a < 2)) when plotted in a log- log scale, namely the emergence of an inflexion point 
before the straight line segment [29]. We must emphasise that this similarity by no means 
implies a possible application of Levy-Gnedenko generalisation of the CLT which states 
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Table 1: Values of the maximum value of P PCzv) maX ) and r] parameters characterising 
PDF Eq. (jHl) for positive, rj + , and negative, rj-, branches for each value of N calculated. 



P{X N ) 



max 



10 



0.350 
0.203 
0.141 
0.128 
0.118 
0.109 
0.101 



10 2 

10 3 
2 x 10 3 
4 x 10 3 
8 x 10 3 



1.80 
1.57 
1.52 
1.47 
1.41 
1.30 



1.81 
1.57 
1.52 
1.47 
1.41 
1.29 



1.6 x 10 4 



that, the probability density function of the sum of ./V variables, each one associated with 
the same distribution Eq. (jHJ) (0 < r\ < 2 and u G 3?) converges in the limit as N goes to 
infinity to a a-stable Levy distribution with a = rj. It is easy to verify that the results we 
report in this manuscript do not follow this theorem because; i) x t and y t variables are 
boxed up within intervals from —1 to +1 and hence their standard deviations are always 
finite, ii) due to the weak chaotic properties, variables x t and y t are not idenpendent at 
all instants t. 

Trying to infer about the scaling behaviour of P (^"jv) max with N, a clear-cut power-law 
behaviour could not be found as it is visible in Fig. [5J 

In the absence of clear power-law behaviour and using the fact that rj decreases as N 
increases, we have tried to extrapolate a value of r\ (N — > oo). To this end, inspired by 
finite-size scaling relations of critical behaviour [27], we have used the following ansatz, 



From this, we have obtained rj^ = 1.02±0.06 (see Fig. H]) very close (within error margins) 
to the exponent of the Lorentzian distribution, Ci (X N ). 

4 Final remarks 

In this manuscript we have presented a numerical experiment on the addition of de- 
terministic variables generated by a conservative map, the triangle map of Casati and 
Prosen [16] . The study has been performed in two different regimes, case I and case II, by 
iterating the map from a set of initial conditions which have uniformly been placed within 
interval [—1, 1). In case I, for which the map is ergodic and mixing, the outcoming stable 
PDF is the Normal distribution for both and Yjv, in perfect accordance with standard 
theory. In case II, for which the map is weakly ergodic for sure and with apparently no 
mixing, we have obtained PDFs which are well described by power-laws for large values 
of the variable. Moreover, the parameter characterising the PDF is smaller than 3 and it 
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presents a decreasing trend as the number of variables N are augmented □. Our results 
are to some extent surprising seeing that, notwithstanding the map is weakly ergodic, it 
fills the phase space as time evolves. Accordingly, it would be expected that Gaussian 
behaviour would be approached (as in case I) as iV increases and not the opposite as 
we have reported here. In such a scenario of weak chaoticity, the map was expected to 
present some anomalous (quasi-steady) behaviour before total occupancy of the phase 
space took place which would imply a crossover to the Normal distribution. This last 
description is analogous to what has been observed in the fractal dimension, which tends 
to increase towards the Euclidean dimension, for the case of low-dimensional sympletic 
maps (see details in Ref. [28]). However, for the case we have studied, the observed in- 
creasing departure from the Gaussian distribution with increasing iV points to another 
direction^. Last of all, from the ansatz [TTJ we have extrapolated the value of r\ when 
N — > oo, which has shown to be 7700 ~ 1, i.e., the same decay as the Lorentz distribu- 
tion, which corresponds to a g-Gaussian with q = 2 in the non-extensive formalism. In 
defiance of the different nature of both systems, our result for the Casati-Prosen conser- 
vative system has provided a similar qualitative result as obtained in Refs. [T3] for the 
logistic (dissipative) map. In other words, departing from variables with a finite standard 
deviation that evolve according with a dynamical system in which strongly ergodic and 
strongly mixing features are not verified, we have been able to define a new variable whose 
limit distribution (spanning the whole domain of real numbers) has a different attractor 
that for the Gaussian. Moreover, in our results the mixing property looks to be a crucial 
element. It is our expectation that, along with results for dissipative dynamical systems, 
this work should suggest further studies on the introduction of an analytical framework 
for cases where the Bernouilli shift is not verified. 
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Appendix 

The method introduced by Meerschaert and Scheflier [25] is based on the asymptotic limit 
of the sum of the variables of dataset {Ajv} under scrutiny. For heavy tail data these 
asymptotics depend only on the tail index of the probability density function, and not on 

3 An exponent equal to 3 corresponds to the lower bound for finite second order moment of distributions 
defined between —00 and 00. 

4 TV can also work as a measure of time. 
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the exact form of the distribution. Hence, if X elements of a dataset are identically and 
(in)dependently distributed, and in addition its probability density function presents an 
asymptotic behaviour, 

P(X N ) ~ \X N \- V ~\ (\X N \^oo), 
it can be proved (Theorem 1 in Ref. [25]) that, 

1 _ ln + [S? =1 (X N , - (X)) 2 ] 

n ~ 21nJ ' (A1) 

where (X) is the simple average (X) = X~ 1 Yj L i=1 X N , i and ln + [x] = max {In x, 0}. 
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Figure 1: Upper panel: Probability density function P (X' N ) vs X' N for scaled variables 
X' N = Xn — (Xn) /cjv, obtained in case I where (Xn) represents the average of Xn [on 
log-linear scale]. Lower panel: The same as the upper panel, but for variable Yn- In both 
panels the line labelled Gaussian corresponds to Eq. (jHJ) with <tn = 1. 
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1000 



Figure 2: Standard deviation, a N , of X N vs. N for case I [on log-log scale]. The fitted 
straight line has a slope of 0.49 ± 0.01. 
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Figure 3: Standard deviation vs N of obtained from map (CQ) [on log-log scale]. 
The dotted line corresponds to Eq. [7] and the grey dashed line has slope 0.94 ± 0.02 
corresponding to the best linear fit on this particular scale. 
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Figure 4: Values of the exponent for positive side, rj + , vs. iV 1 [on log-log scale]. The line 
represents a numerical adjustment for Eq. ffTTj) with r/^ = 1.02, c = 0.049, and 5 = 0.40 
( x 2 = 6.9 x 10" 4 and R 2 = 0.986). 
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Figure 5: Maximum value P pGv) max vs. iV 1 according to the values of Table Q] [on 
log- log scale]. 
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Figure 6: Probability density function P' (X N ) = P ^x N ) vs [in log-log scale]. The 
number of initial conditions (points used to construct PDFs) is 10 7 except for iV = 8000 
(7.5 x 10 6 points) and iV = 16000 (7 x 10 5 points). The power-law decay is evident for 
large X^. The two straight lines have slopes — (r] + + 1) with rj the exponent for iV = 1000 
and N = 16000 (see Tabled]). 
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